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Critical surface in hot and dense QCD with the vector interaction 
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We discuss tlie chiral phase transition of hot and dense quark matter. We illustrate that the first- 
order phase transition is generally favored at high baryon density and the repulsive vector-vector 
interaction weakens the first-order phase transition. We use the Nambu-Jona-Lasinio model with 
the Polyakov loop coupling for concreteness. We locate the QCD critical surface on the quark mass 
plane for various values of the vector coupling constant. We find that, with increasing quark chemical 
potential, the first-order region in the quark mass plane could shrink for sufficiently large vector 
coupling. This may be a possible explanation for the recent lattice QCD results by de Forcrand and 
Philipsen. 

PACS numbers: 12.38.Aw, ll.lO.Wx, ll.30.Rd, 12.38.Gc 



I. INTRODUCTION 

On the phase diagram in the plane with the axes by 
temperature, density, pressure, concentration, external 
fields, etc, the "critical point" commonly refers to the 
terminal point of the phase transition boundary of first 
order. This point has an exact second-order phase tran- 
sition regardless of symmetry properties. 

Interestingly enough, the existence of the critical point 
is a possibility not only in condensed matter physics but 
also in QCD (Quantum Chromodynamics) physics 
The QCD phase diagram with the temperature T and 
the quark chemical potential has taken on significance 
in the application to the heavy-ion collisions at various 
energies as well as to the neutron star structure. There 
seems to be a consensus that the QCD phase transitions 
associated with chiral symmetry restoration and color de- 
confinement are smooth crossover when T goes up with 

~ 0. In contrast, the situation changes in the different 
regime where ^ grows with T ~ 0. Model studies 0, Q 
suggest that the first-order chiral phase transition should 
occur at some ^ comparable to the constituent quark 
mass (or one third of the baryon mass) when T ~ 0. To 
be consistent with crossover at T 7^ and fi — 0, there- 
fore, the critical point is expected to exist at intermediate 
Te and /i^. 

It is, however, important to remark that the existence 
of the QCD critical point is still under dispute [J. One 
interesting negative observation against the QCD critical 
point has come from the Monte-Carlo simulation of QCD 
on the lattice with varying the light-quark mass m„d and 
the strange-quark mass 5], though some other lat- 
tice simulations are rather affirmative. (For a review see 
Ref. 6]). 

For the symmetry reason the chiral phase transition 
in three-flavor quark matter at m„d = i^s = is pre- 
sumably first order [7] . Because the mass term explicitly 
breaks chiral symmetry, the first-order transition turns 
to crossover at some m„d and Q, which defines 

the critical boundary in the m„d-ms plane. Model stud- 
ies that support the QCD critical point predict that the 
first-order region in the TO„d-ms plane expands by the 



effect of increasing ji, so that the physical quar k mass 
point hits the critical surface at ^ = llll|. How- 

ever, de Forcrand and Philipsen 0] recently claim that 
the first-order region should not expand but rather shrink 
at higher density as long as ^/T is small. 
The purpose of this paper is twofold: 

1) We will extract the general feature of the quasi- 
particle description at high density to favor the first-order 
phase transition. We make use of simple closed expres- 
sions by limiting our discussions to the T = case. In 
fact, the confirmation of the first-order phase transition 
at r = and /i 7^ suffices for the existence of the QCD 
critical point since the QCD phase transition at /i = 
and T ^ is crossover. In the same way we discuss the 
effect of the repulsive vector-vector interaction on the 
general ground. These are all discussed in Sec. [TTl 

2) We will point out in Sec. IIIII that the shrinkage of 
the first-order region in the rriud-fTis plane is not uncom- 
mon once we take account of a vector- vector interaction. 
We draw the critical surface using the Polyakov loop aug- 
mented Nambu-Jona-Lasinio (PNJL) model to exemplify 
the effect of the vector interaction, which may be a likely 
explanation for the results by de Forcrand and Philipsen. 



II. FIRST-ORDER PHASE TRANSITION 

Let us first consider a simple chiral model at T = 
with and without the vector interaction. We will do so 
because we should understand the underlying mechanism 
for possibility of the first-order phase transition at T = 
and / to elucidate the effect of the vector interaction. 

The simplest case in the chiral limit is enough for our 
present demonstration in which one dynamical mass M 
serves as the order parameter. Here M is either the con- 
stituent quark mass in the color deconfined phase or one 
third of the nucleon mass minus binding energy of nu- 
clear matter if the system confines color. Although we 
can formulate both, we shall focus on the former, i.e. de- 
confined quark matter, hereafter. (For the latter, see dis- 
cussions in Ref. ,3j.) We assume that the pressure ^'^[M] 
results in the chiral symmetry broken phase at zero den- 
sity. That is, the free energy, — P^[Af], has minima lo- 
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cated at M 
form; 



±Afo; we simply postulate the following 



P^[M]^-a{Ml-My 



(1) 



with a parameter a. Here we note that a linear term in M 
should be present if the current quark mass is nonzero. 
We can neglect this explicit chiral symmetry breaking in 
the qualitative level because such a term has only minor 
effects on the phase transition in the two-flavor sector. 
In the three-flavor case, in contrast, the Ua(1) break- 
ing term generates a cubic term in M which favors the 
first-order phase transition. We will not think of this 
situation; our purpose here is to see how the first-order 
transition is possible at high density even though it is of 
second order at vanishing density. Thus, the above form 
of Eq. ^ is valid when all the quarks are massless and 
the three-fiavor Ua(1) breaking is not significant. 

Now let us turn finite on. As long as /i is smaller 
than the lowest-lying mass of fermionic excitation, noth- 
ing happens and the vacuum remains empty. Once /i 
exceeds a certain threshold M, a finite amount of den- 
sity appears. The pressure has a contribution from the 
density which is generally expressed as 



dfi'ntj{fj.') . 



(2) 



Here nq{^) represents the fermion density. In the quasi- 
particle picture it is given by the integrated Dirac-Fermi 
distribution function with the constituent mass; 
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(3) 



where e = yj5^ + M^ and v is the fermionic degrees of 
freedom (color x flavor x spin). In two-flavor quark mat- 
ter, for relevant example, z/ = (3 colors) x (2 flavors) x 
(2 spins) = 12. We note that 6 denotes the Heaviside 
theta function, which signifies that the system at n < M 
is empty. In fact, the theta function is essential to make 
a double-peak shape in the total pressure, as we will see 
soon. 

It is possible to perform the integration ^ to find an 
analytical expression with logarithmic terms. To simplify 
our qualitative analysis, however, we shall introduce an 
approximation as 



247r2/i2 



(4) 



which turns out to be a good approximation as shown 
in Fig. [T] The solid curve represents Eq. ([l]), while the 
dotted curve is Eq. ^ with Eq. ([3]) substituted. Because 
more particles can reside in the Fermi sphere for smaller 
mass, Pi^i\M] has a maximum at M = 0. 

Let us consider the condition for P\M] = P^[M] 
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FIG. 1: Comparison between the exact integration in Eq. ((2)) 
(by the dotted curve) and the approximation in Eq. (by 
the solid curve). 
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FIG. 2: Sketch of the double-peak pressure P[M] resulting 
from the sum of Px[M] and Pf,[M]. 



and have a peak at M = Mq and A/ = re- 

spectively (see Fig. [5]). The existence of double peaks 
in P[M] requires that ^ < Mq, meaning that /i should 
not be much greater than Mq. [So, fi can be larger than 
Mq slightly.] This is necessary for the peak at M = Mq 
to survive. At M = the pressure curvature (i.e. the 
coefficient of the term) should be negative, that is; 



a < 



167r2 M§ 



< 



167r2 



(5) 



At the first-order critical point the peak at Af = is as 
high as the second peak at M = Mq (neglecting a small 
shift by the contribution from P^[M]), which yields the 
critical condition that 



247r2 M^ 



(6) 



As long as ^ is raised with fi < Mq satisfied, the curvature 
condition ([5]) is sufficient for the existence of /ic deduced 
from Eq. ([6]). This is another way to see why we should 
have required fi < Mq. 

We shall next take account of the mean-flelds from 
P^[M] to have a flrst-order phase transition. Here P^[M] the vector- vector interaction, — G'v'('07^'(/')('07''V') [Hlj 
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[H, [T3| ■ Because {'ip^^il') is nothing but the fermionic 
density, roughly speaking, the vector interaction gener- 
ates a contribution to the pressure; 



Pv[M] 



-Gynl 



M'^f e{^i^ - M'"-) (7) 



at T = 0, which takes the same functional form as the 
approximated (|3]). [In the above we have constructed 
Pv[M] using Uq given by dPfj^/djjL. This approximation is 
qualitatively reasonable, but not self-consistent with the 
full pressure from which nq should have been inferred. 
We will come back to this point later.] The coefficient in 
Eq. ([U is hence modified by the effect of Eq. ([7]) (that is, 
the effective fermionic degrees of freedom are reduced), 
and the curvature condition is then 



a < 1 



37r2 



167r2 Af 2 



(8) 



Therefore, even though we start with a that satisfies 
Eq. (I5|), there is a critical Gv for which Eq. ([8|) is not 
satisfied within /i < A/q. Then, the first-order phase 
transition would disappear. 

Let us see concrete numbers; in the conventional NJL 
model with two flavors [l^, for which our discussions 
based on Eq. ([T]) are valid, one can read a as 
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0.067, 



(9) 



where we used v = 12 and the ultraviolet cutoff A = 
631 MeV, the four-fermion coupling GgA^ = 2.19, and 
the resultant Mq = 336.2 MeV [H]. The first term orig- 
inates from the zero-point energy and the second from 
the four-fermion interaction. Then, the right-hand side 
of Eq. ([5]) is 0.076 and the inequality is certainly satisfied. 
In the case of the chiral quark model (l6| . for compar- 
ison, a is estimated as a ~ m'^f'^/ {8Mq), which yields 
a = 0.02 ~ 0.05 depending on the value of the a meson 
mass. This is again within the region of the first-order 
phase transition. 

The critical point in the NJL model case is 



= 1.07Af, 



: 



(10) 



from Eq. ^ . This estimate is consistent with the empir- 
ical value in the NJL model study. 

It is easy to see the effect of Gv from Eq. ([5]). As we 
have confirmed the critical /i is nearly Mq and we can 
replace /i Mq in Eq. ([S]) approximately. Then, the 
inequality does not hold when 



Gv > 0.25Gc 



(11) 



meaning that such Gv makes the first-order phase tran- 
sition disappear at all along the /i axis. It is impressing 
that this rough estimate is fairly consistent with a more 
serious analysis in the NJL model [l^l . 

We can learn from the above argument that the repul- 
sive vector- vector interaction reduces the pressure arising 
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FIG. 3: Critical surface with standard curvature in the PNJL 
model without the vector interaction. The physical mass 
point in this model is {niud ~ 5.5 MeV, tUs = 135.7 MeV) 
which hits the critical surface at /i = he- 



from the degenerated particles in the Fermi sphere. To 
achieve chiral restoration with such strong vector inter- 
action, then, /i > Mq is necessary and a peak around 
A/ = Mo is washed away by the tail in P^[Af] -I- Pv[M] 
which extends up to Af = /i. This is the qualitative 
mechanism how the vector interaction would hinder the 
first-order phase transition. 



III. THREE-FLAVOR MASS PLANE 

From the point of view of thermodynamic pressure, 
so far, we have seen how the first-order phase transition 
could occur at ^ = — especially in the case of 
two-flavor cold quark matter. 

Let us now proceed to the main part of our discus- 
sions on three- flavor quark matter. For a fixed value of 
/i, the chiral phase transition at finite T is of first order 
when the current quark masses rriu^ and stay small, 
and becomes of crossover if m^^; and mg exceed a crit- 
ical boundary. Therefore, in the 3D space of mud, i^s, 
and /i, the critical condition makes a hypersurface. For 
concreteness we adopt the PNJL model used in Ref. [ll| 
to draw the critical surface. The NJL model might al- 
low for artificial quark excitations at finite temperature 
due to the lack of confinement. This tends to shift the 
location of the critical point {ij,e,Te) to higher he and 
lower Te because the temperature effect smears the first- 
order phase transition induced by density. The PNJL 
model cures this (for instance Te changes ^ 50 MeV in 
the NJL model to ~ 100 MeV in the PNJL model with 
the same parameter set [l3|) by means of the color pro- 
jection by the Polyakov loop which is a colored phase 
factor associated with single-quark excitation. Since the 
model details are not our main focus, we simply refer to 
Ref. [ll[i3 for details. 

Figure [3] is the critical surface with standard curvature 
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obtained in the PNJL model. Above the critical sur- 
face the finite-T phase transition is of first order, while 
no sharp phase transition takes place below the surface. 
The second-order phase transition sits on the critical sur- 
face. This is just one model example but shows typical 
behavior in model studies The physical mass point 
hits the critical surface at = = 313.5 MeV where 
the critical temperature is T = Te = 101.8 MeV (shown 
by a cross in Fig. [3|). 

It is general in the quasi-particle picture that the den- 
sity effect induces a pressure like Eq. ^ whose maximum 
is located in the (partially) chiral symmetric phase. In 
the presence of such an additional peak near Mud — nT-ud, 
Ms ~ TTis, the first-order transition region is enhanced 
with increasing /i. 

The problem is that this curvature of the critical sur- 
face as shown in Fig.[3]might be inconsistent with the lat- 
tice observation [Sj] even though the model results seem to 
be rather robust not relying on any special assumption. 

We here propose one scenario that has a natural ac- 
count for the lattice results and, at the same time, 
may not necessarily exclude the existence of the criti- 
cal point. The necessary ingredient is the vector- vector 
interaction alone. As a matter of fact, the repulsive 
vector interaction is anticipated from the hadron prop- 
erty [13, [3 • Figures S] and \E\ arc the examples from the 
PNJL model with the vector interaction incorporated. 
Here we should explain the mean-field treatment for the 
vector interaction in a self-consistent way 0, [ll|, [l^, [H, 
The zeroth component in the vector interac- 
tion, — Gy(Vj7o'!/')(V'7°V'); produces the mean-field terms, 
-2Gyn,V57> + Gvnl 
chemical potential as fi 
justs the larger contribution from the former. Then, the 
mean-field nq is fixed by dP/duq — which guarantees 
the thermodynamic relation; nq — dP/dfi = 0. [Note that 
dP/duq = -2GvdP/dti.] 

We have chosen a considerably large value of Gy to 
make it easier to perceive the effect visually. The point 
is that the inclusion of Gy induces the opposite 
curvature to that in Fig. [3] at small fi, and eventually 
the curvature returns to the standard one at larger ^. It 
should be noted that the first-order region at ^ = does 
not change but the axis scale in Figs. [Hand [S] is different 
from Fig. O 

Although our aim is just to present some demonstra- 
tions like Figs. [4] and [5l it is intriguing to make a quan- 
titative comparison. Along the diagonal m = rriud = "m-s 
line, the critical mass is expanded in terms of fi as 



The former term adds to the 
■ fi — 2Gvnq and the latter ad- 



m,(0) 



= 1 



C4 



JJ_ 



(12) 



Our calculations give C2 — 5.88 and C4 = 43.8 at Gy = 
as shown by the solid curve in Fig. [HI The lattice results 
are, on the other hand, C2 = —0.7(4) and C4 ~ in the 
first paper of Ref. 0] which is indicated by the upper 
shaded region in Fig. [51 In the last paper of Ref. [Hi sig- 
nificantly different values are reported; C2 = —3.3(3) and 
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FIG. 4: Demonstration of the effect of the vector interaction 
in the PNJL model with Gv = OAGs- 
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FIG. 5: Demonstration of the effect of the vector interaction 
in the PNJL model with Gv = O.SGs- 



C4 = —47(20) as shown by the lower shaded region. In 
our calculations the vector interaction drastically alters 
the curvature from the solid curve to the dashed (dotted) 
curve in Fig.[n|at Gy = O.4G5 (Gy = O.8G5). 

From Fig. [H] we may say that some Gy < O.4G5 could 
be enough to understand the results in the first paper of 
Ref. , while the vector interaction alone is not sufficient 
to reproduce the results in the last paper of Ref. . To 
make the statement more conclusive we need more lattice 
QCD data. 

The rest of this paper will be devoted to qualitative 
explanation of this back-bending curvature as a result of 
the vector interaction. We will discuss the mechanism in 
order. 

1) We have to realize that the finite T is important to 
understand the back-bending behavior. The information 
of the critical temperature is implicit in Figs. [3l [H and 
[5] The general tendency is that the critical temperature 
becomes smaller as increases. Thus, in fact, the critical 
surface is cut at large chemical potential for which the 
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FIG. 6: Critical mass nic along the rriud ~ ms line as a func- 
tion of fj,. The (red) solid curve represents the PNJL model 
results at Gv ~ 0, the (green) dashed curve at Gv = OAGs, 
and the (blue) dotted curve at Gv = O.SGs- The upper and 
lower shaded regions correspond to the lattice results reported 
in the first and last papers in Ref. respectively. 
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FIG. 7: nl (solid curve) and (dotted curve) at T/fi — 
and T/fi — 0.3, where is normalized such that and 
nearly coincide at T/jj. = 0. 



critical point touches T = 0. 

2) In the vicinity of /i = 0, hence, the finite T effect 
is substantial, which makes the functional forms of the 
integrated nq{^) with respect to /i (i.e. P^) and de- 
viate from each other unlike the T = case. It should 
be mentioned that the simple argument in Sec.|lT]is stiU 
apphcable for m„d — rus as it is (except for the Ua(1) 
breaking term) . Then the integrated one ^ turns out to 
be less steeper than as a function of M, which can be 
readily confirmed by simple numerical calculations (see 
Fig. ED. 

Thus the pressure contribution by density-induced par- 
ticles is relatively more suppressed by finite T than the 
vector interaction effect involving n^. This means, in 
other words, that brings a sharper modification than 
to the total pressure at higher T (and thus smaller ji 
in turn). 

3) Because the repulsive vector interaction has the op- 



posite effect to the density, as we have seen in Sec. |TT1 
some Gv exists which is large enough to invert the den- 
sity effect. Still, the first-order region at /i = is intact 
since = at /i = 0. 

4) With appropriate Gy 7^ 0, as ^ goes up, the first- 
order region shrinks by the effect of growing from the 
vector interaction which overwhelms the effect of as 
long as ^./T is small. 

5) When ^ gets larger and T//i becomes smaller, the 
functional shape of comes to be identical to that of 
P^ as seen in Fig. [71 If Gv is not too large, the den- 
sity effect can surpass the vector interaction eventually, 
which makes the curvature bend back into the standard 
direction at high which is manifestly the case in Figs. [3] 
and[5j 



IV. SUMMARY 

We have clarified how the quasi-particle description 
can lead to a first-order phase transition in cold and dense 
quark matter. In the same way we have intuitively made 
clear the role played by the mean-fields from the repul- 
sive vector- vector interaction which reduces a pressure 
contribution from density-induced particles. 

We have then discussed the vector interaction in the 
mud-fTT-s plane with three fiavors. We have drawn the 
critical surface using the PNJL model with the vector 
interaction. Along the same line as the quasi-particle 
description we have given a simple account for the back- 
bending behavior of the critical surface, which is con- 
sistent with the negative curvature in the recent lattice 
simulations. Logically speaking, therefore, there might 
be a chance that the critical point still persists even with 
the negative curvature at small density. 

It should be remarked, however, that whether the back 
bending occurs or not cannot be naively interpreted as 
whether the critical point exists or not. This is because, 
as we mentioned, the critical surface is cut at some large 
/Lt where the critical temperature is zero. As a matter 
of fact, the critical value of Gy is around 0.25 as in 
Eq. pT|) (see also Refs. [O, [H, [13) for which we cannot 
observe clear back-bending behavior in the PNJL model. 
In this way, if we take the comparison to the PNJL model 
analysis seriously, the lattice results by de Forcrand and 
Philipsen certainly suggest the nonexistence of the QCD 
critical point. 

There are important issues to be investigated in the 
future. First of all, the determination of Gy is indispens- 
able. For this purpose, unfortunately, the description of 
the vector meson property within the NJL model is not 
quite reliable because the momentum cutoff is not large 
enough as compared to the vector meson mass. The lat- 
tice simulation is one possibility, though observable to 
measure Gv is not clear yet. Second, not only the vector 
interaction but also the anomaly ('t Hooft) term is impor- 
tant for the existence of the QCD critical point Q. The 
density dependence of the anomaly coupling strength is 
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unknown from the lattice simulation, which may affect 
the critical surface curvature. The density dependence 
of Gv is not known, either, but the point of what we 
have shown here is that intriguing density dependence 
appears from even for a constant Gy ■ Hence, if the 
opposite curvature of the critical surface to the standard 
one is really established in the lattice simulation, it could 
be interpreted as a circumstantial evidence for finite re- 
pulsive vector-vector interaction. 
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